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Abstract 

A method is developed to consistently satisfy the Gibbs equilibrium conditions 
between the quark-gluon and hadronic phase, although each phase has been formulated 
in separate grand canonical partition functions containing three quark flavours. The 
sector in the space of thermodynamic variables, where the transition takes place, is 
restricted to a curve, according to the phase diagram of QCD. The conservation laws of 
quantum numbers are also imposed on the transition curve. The effect of the inclusion 
of the newly discovered pentaquark states is considered. The freeze-out conditions of 
S + S, S + Ag (SPS) and Au + Au (RHIC) are found compatible with a primordial 
QGP phase, but the conditions indicated by Pb + Pb (SPS) are not. 
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1. Introduction 

Quantum Chromodynamics is universally accepted as the theory of strong interactions. 
Within the context of this theory the phase of Quark- Gluon plasma receives accurate de- 
scription. However, the formation of the hadronic phase, which is the final state of any 
possible primordial QGP state, still remains an open problem in view of QCD. On the other 
hand, the hadronic multiplicities emerging from heavy-ion collisions have been extensively 
and successfully predicted by statistical models using a handful of thermodynamical param- 
eters [1-8]. So the use of two separate models for the QGP and the hadron phase, called 
Hadron Gas (HG), offers a complementary approach. 

QCD predicts that the transition between QGP and the hadronic phase is a first order 
one at high baryon densities (depicted by a critical line on the (T, yU^) plane), while it is of 
a higher order at small or zero baryon densities (crossover). The end point of the first order 
transition line is a critical point [9]. The transition points must be restricted to a curve on 
the phase diagram of temperature and baryon chemical potential. In view of this aspect any 
models to be used for the description of QGP and HG have to be matched properly at the 
transition between the two phases. 

The aim of this work is to trace the sector of the space of thermodynamic variables, where 
the QGP-hadron transition occurs, with the following requirements: a) Any mixed phase 
formed in the first order part of the transition must occupy only a curve in the space of the 
thermodynamic variables. This requirement is even stronger in the crossover area where a 
mixed phase does not exist, b) The Gibbs equilibrium conditions have to be satisfied, which 
amount to Tqqp = Tjjg for thermal equilibrium, Pqgp = Phg for mechanical equilibrium 
and {ij>}qgp — {a*}hg for chemical equilibrium, where {/i} stands for the set of chemical 
potentials used in the description of the two phases, c) All the conservation laws of quantum 
numbers hke baryon number B, electric charge Q, strangeness -S", etc. have to be satisfied 
at every point on the transition line, in a way that they could be extended for every number 
of flavours that are present in the system. 

These problems are confronted every time separate partition functions are used for the 
two phases, but the simultaneous fulfilment of the above conditions has not been achieved. 
Among the numerous examples that exist, in [10], where only light, identical quarks are used 
{u = d = q), the curves of equal pressures arc made to approximately coincide by a choice 
on the external parameters B and a^, something which does not allow matching when other 
flavours are introduced. In [11] the strange fugacity is discontinuous at the HG-QGP 
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transition and the conservation of baryon number can only be accommodated in the case 
of first order transition. In [12] the strange chemical potential fig is also discontinuous. In 
[13] only q quarks are considered and the requirement of continuity of chemical potentials 
and conservation of the baryon number leads to a mixed phase which occupies a surface 
and not a line on the {T,iib) plane. The same is true in [14-16] where also s quarks are 
included. In [17] there is an analogous situation as in [13] with a critical line at the (T, /i^) 
plane but the conservation of baryon number is not considered. In [18] the q and s quark 
chemical potentials are continuous but baryon number and strangeness of the system are not 
kept constant during hadronisation, since hadrons evaporate from QGP. The considerations 
of [11-18] are consistent with a first order transition but cannot be valid at the crossover 
region. In this work all the thermodynamic variables and the pressure will be kept continuous 
at the transition line (in contrast with [10-12]), the first order part of the transition will be 
presented by a line on the {T,iib) plane (differing from [13-16]), in the mixed phase the 
quantum numbers will be conserved to each constituent phase (differing from [14]) and no 
evaporation of hadrons will be assumed from the system (differing from [18]). 

Let us consider the requirements that a system with Nf quark flavours has to satisfy. 
Every conservation law accounts for two equations to be fulfilled. One sets the value of the 
quantum quantity, e.g. {B)qqp — hi and the other assures the conservation at the phase 
transition, e.g. {B)qqp — {B)jjq. The total number of equations that must hold are, thus, 
2Nf + 1 (the unit accounts for the equahty of pressures). Assuming the existence of full 
chemical equilibrium, every quark flavour introduces one extra fugacity in the set of the 
thermodynamical variables, which, with the inclusion of volume and temperature, amount 
to Nf + 2. At the crossover region, the surviving free parameters required to fulfil the 
necessary equations decrease to Nf + 1, because of the equality of densities and consequently 
the equality of volumes between the two phases {Vqgp — Vhg)- At the first order transition 
line the free parameters are Nf + 2, since now Vqgp 7^ Vhg- It is evident then that the 
necessary 2Nf + 1 conditions can be fulfilled only at the first-order part of the transition and 
only when there is one flavour present, Nf — 1, ot when the u and d quarks are considered 
identical {q quarks, described by a single chemical potential fj,q). It has to be clarifled that 
the conditions like {B)qqp — hi have to be satisfled in order to have a whole line of transition 
points. If these equations arc dropped, then we arc left with Nf + 1 equations, which can be 
solved, but result to a unique point in the space of the thermodynamical parameters. 
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2. Expanding the fugacity sector 

It is clear that in order to satisfy 2Nf + 1 relations, every flavour has to be accompanied 
by two fugacities instead of one. The multiplicity data emerging from heavy ion collisions 
suggest that the thermalised hadronic system has not achieved full chemical equilibrium. 
First the strangeness partial chemical equilibrium factor 7^ had been introduced [2] and 
used extensively to model the data [3-4]. Also a similar factor for the light quarks 7^ was 
introduced [5] and used in certain analyses [6]. Here the light u and d quarks will be 
accompanied by separate fugacities 7^, 7^. A factor 7^- controls the quark density rij+nj 
in contrast with the usual fugacity \j which controls the net quark density rij — rij [3] . 
These additional fugacities can serve the purpose of satisfying the necessary equations at 
the transition point, as well as, preserving the continuity of chemical potentials between the 
two phases. 

A system with 3 flavours {u, d and s quarks) is described by the set of thermodynamical 
variables (T, A^, 7„, A^, jd, K, 7s) = {T, {A, 7}). Setting x = Vhg/Vqgp, the set of equations 



to be satisfied at every phase transition point will be 

^QGP 

(r,{A,7})=P,,(T,{A,7}) (1) 

UB^apiT, {A, 7}) = X UB^aiT, {A, 7}) (2) 

i^Qqgp (T, {A, 7}) = a; ng^^ (T, {A, 7}) (3) 

ns^apiT, {A, 7}) = X ns.aiT, {A, 7}) (4) 

{T, {A, 7}) = 2/3 uqq^p (T, {A, 7}) (5) 

n5^^,(r,{A,7}) = 0, (6) 
where n denotes densities. For isospin symmetric systems one has to set P — 1 in (5). 



Eqs. (l)-(6) only have one free variable, necessary to produce a whole transition line in the 
phase diagram. At crossover x — 1, whereas at the first order transition line the inequality 
Vqgp 7^ Vhc preserves the survival of x as an extra variable. 

3. A solution for the transition curve 

The above considerations are applicable to every partition function connected to the 
hadronic and the quark state. It is interesting, though, to check whether they can produce 
a real solution for the transition curve, i.e. that the system of equations (l)-(6) is not im- 
possible. For this reason two simple models, will be employed. For the Hadron Gas phase 



only the repulsive part of the hadronic interaction will be taken into account through a Van 
dcr Waals treatment of the system volume. Also the correct Bose-Einstcin or Fcrmi-Dirac 
statistics applicable to each hadron will be considered. The hadronic partition function will 
be extended to all hadronic states containing u, d and s quarks as they are listed in [19]. 
First, the HG partition function for point particles can be written down as 

Y /"oo 1 

lnZHGpt{V,T,{X,-f})^—^Y.T.9.^ I ^ JM/T, , ' 

where gg,i are degeneracy factors due to spin and isospin and a — —1(1) for bosons (fermions). 
Index a runs over all hadronic families, each of which contains members with the same quark 
content and i over all the particles of this family. The fugacity Aa — Ylj X^^ ^j^^^+^j ^ where 
j = u,d,s and Nj{Nj) is the number of j{j) quarks contained in a hadron belonging to family 
a. For the light unfiavoured mesons with quark content {ci/2){uu + dd) + C2ss, Ci + C2 = 1, 
the fugacity used is Aa = {luldY^ll"^ ■ 

If each hadron i is assumed to occupy volume T^, then the available volume for the system 
is reduced to IS. — V — Y^^NiVi, where Ni is the number of particles i present in the system. 
Assuming that each hadrons' volume is proportional to its mass, then Vi/rrii = Vq, with 
Vo remaining an open parameter controlling the hadron size. The available volume can be 
written ets A — V — J2i NimiVo. Dividing by A we get 

V . V 
1 = X ~ Ppt^o 



A '^^^ " 1 + pptVo ' 

where Ppt is the mass density of a system of point particles at volume A. Defining the 

quantity / = PptVo, it can be calculated to be 

where the sum over all particles is organised first to a sum over all families. The available 
volume is then 

The partition of the extended hadrons at volume V may be considered equal to the 
partition function of the point particles at volume A 



The pressure is 



Pi 



In Z{V, T, {A, 7}) = In Zpt{A, T, {A, 7}) ■ (10) 
dlnZHciV,...) din ZhgAA,...) 



HG 



dV dV 

A 



or by using eq. (9) 



PHa = ^. (11) 



In a similar manner the densities are calculated to be 

UiHG^ (12) 

For the QGP phase a simple model containing 3 flavours is used. The quarks are non- 
interacting and only the presence of gluons is accounted for, as well as the effect of the vacuum 
through the MIT bag constant, B. A wealth of quark fugacities is easily accommodated, 
though, in this model. The QGP partition function is consequently 

N N V r°° 1 Stt^T^ BV 

(13) 



where Ng = 2 and = 3. Index j runs to all quarks and antiquarks and the fugacity 
Xj — Xj^ and 7j = 7^-. The current quark masses are m„ = 1.5, = 6.75 and = 117.5 
MeV [19]. 

At the first order QGP-HG transition a mixed phase is assumed. This phase spans over 
a curve in the (T, {A, 7})-space, so these variables are kept constant throughout the mixed 
phase. The only thermodynamic variable which is allowed to change is the system volume V. 
The volume equals Vhg at the pure hadronic phase, at one end of the first order transition 
and Vqop at the pure quark phase, at the other end of the transition. The partition function 
of the mixed phase can then be written down as 

In Z^i^ediV, T, {A, 7}) = 5 In Zhg{Vhg, T, {A, 7}) + (1 - 5) In Zqgp{Vqgp, T, {A, 7}) . (14) 

The parameter 5 is < 5 < 1 and ior 5 — 1{5 — 0) we have pure HG(QGP) phase. Since 
in all cases the volume appears as a multiplicative factor in the partition function, the ratio 
g = does not depend on the volume. Then eq. (1) suggests that at every transition point 

QHGiT, {A, 7}) = 9QGP{T, {A, 7}) . (15) 

Then eq. (14) is equivalent to 

In Z^i,ed{V, T, {A, 7}) = [6Vhg + (1 - S)Vqgp] 9hg{T, {A, 7}) , (16) 

suggesting that the system volume V of the mixed phase can be defined as 

V^5Vhg + {1-S)Vqgp. (17) 



Through eqs. (1), (16) and (17), the pressure in the mixed phase can be calculated for every 
V to be 

p _ kT\uZmixed{V,T,{X,-i}) (rs\^\\-p -p 

^mixed — y — Kl gHG[-L , {A, 7|j — -Tf/G — ^QGP ■ (.J-oj 

The pressure is, consequently, kept constant throughout the first order transition. The part 
of the P — V isotherm which corresponds to the mixed phase is parallel to the V axis, as it 
is expected in a first order transition. 

On the contrary, the densities vary. Let Aj be the fugacity associated with the density 
rii. Then 

- J_\ ^^^^niixed{V,T,{X,'y}) _ 

mixed i r i 



1 

V 



V dXi 

c.r X dgHGiT,{X,^}) dgQGpiT,{X,-f}) 
SVuaXi + (1 - ^)yQOP>^i 



or 



5Vhg (1 - 5)Vqgp 

- 5F,o + (l-5)W"^ 5Vno + {l-5)V^Gp''' " ^''^ 
It is easily checked that for 5 = 1(5 = 0) the density of the pure hadronic (quark) state is 

produced. The conservation of the quantum quantities B, Q and S is assured in the mixed 

phase. Eq. (2) implies that {B)jjq — {B)qqp. Then multiplying the two sides of eq. (19) by 

V gives 

VUB mrxed = ^ {B) + (1 " (^)qGP i^) mixed = i^) HG = i^) QGP " (20) 

Eqs. (3)- (4) produce similar equations for Q and S. 

Between the crossover region (where x — 1) and the first order transition line (where 
X 7^ 1) the critical point resides. Observing eqs. (l)-(6) it is clear that they do not pro- 
vide a restriction on x, other than it should be a continuous function. So these equations 

can accommodate an additional constraint in the form of x, which may be provided by a 
sophisticated partition function that records the full part of interaction (attractive part as 
well). 

The system of eqs. (l)-(6) can then be solved for x = 1 for the crossover region or for x 7^ 1 
for the first order transition curve. The system is simplified observing that the strangeness 
neutrality at the QGP phase (eq. (6)) leads to the solution = 1. This solution is valid 
for every case of QGP partition function, as long as products of the fugacities of the strange 
quark with the fugacities of w or 0? quarks do not appear. The initial system then is reduced 



to the system of eqs. (l)-(5) for = 1. The HG partition function (7), (10) and the QGP 
partition function (13) is used to the system of eqs. (l)-(5). For the particular choices of 
partition functions, two parameters, B and Vb, are left open, producing different solutions for 
the transition curves. The system of eqs. (l)-(5) for = 1 accepts as solution for the variable 
75 the value 0, since then eq. (5) is automatically satisfied. This a trivial solution because it 
is equivalent to the absence of the strange quarks in the system and therefore such solutions 
should be excluded. Non-trivial solutions for the thermodynamic variables are depicted for 
the parameters S^/^ = 280 MeV and Vq — 1.4/(4i?) in Figs. 1-4 for the isospin symmetric 
case {(3 — 1). Lines (a) represent the solution without the inclusion of pentaquarks. The 
crossover region, which is determined uniquely after setting the parameters B and Vq, is 
drawn everywhere with slashed lines. 

Two additional matters concerning the phase transition arc the position of the critical 
point where the crossover region ends and the first order transition sets in, as well as, the ratio 
of the volumes x — Vhg/Vqgp at the first order transition line. These two matters cannot be 
dealt with the simple choices of partition functions used for the calculations of this section, 
since the attractive part of the interaction among hadrons or quarks is completely neglected. 
To display certain solutions for the critical curve within the context of the partition functions 
(7), (10) and (13), a position for the critical point has to be chosen. This position is set 
at fjLBcr.p. — 360 MeV, according to the result of [20], where the critical point is located by 
lattice QCD calculations using three quark flavours and considerably reduced light quark 
masses, approaching their physical values. Eqs. (l)-(6) for a; = 1 is a system of 6 equations 
with 7 variables. Adding one more equation (/i^ cr.p. = 360 MeV) a system with 7 equations is 
formed, which is solved to determine the full set of thermodynamic variables that correspond 
to the position of the critical point. 

A form for the ratio of the volumes x also has to be defined. This form has to produce 
a; = 1 at the position of the critical point. Moreover, it is chosen to produce a given value 
Xi at a specific value of A^ i. A simple form which implements these demands is 



where the exponent e regulates the curvature of the first order transition line. For lines (a) 
Xi = 1.1 at A„ 1 = 14.2 and e = 1.5 are chosen. Of course, any function of x can be used, 
producing different first order transition curves. The resulting first order transition lines are 
drawn with solid lines in Figs 1-4, while the respective critical points are represented by solid 
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(21) 



circles. 

Temperature T is displayed as function of the baryon chemical potential /x_b in Fig. 1. 
In the same figure, line (d), which represents the phase transition line as it is calculated 
from the Lattice QCD in [21], is drawn for comparison. The relative chemical equilibrium 
fugacity 7^ is displayed as function of /ib in Fig- 2. This particular solution leads to the 
gradual suppression of ju as baryon chemical potential increases. The connection of 7„ and 
7d, for isospin symmetric solution, is depicted in Fig. 3. The line 7^ = 7^ is also drawn for 
comparison. The relative chemical equilibrium factor 7^ is drawn as a function of the baryon 
chemical potential in Fig. 4. 

In Fig. 5 the ratios of volumes x, which are used in the first order transition, are drawn 
against the baryon chemical potential. The used forms of x are increasing functions with 
respect to the baryon chemical potential. The resulting first order transition lines pro- 
duce smaller temperatures as the baryon chemical potential increases, something which is 
expected. 

One direct consequence of the simultaneous solution of eqs. (l)-(6) is that the relative 
chemical equilibrium fugacities have values that depend on each other at every transition 
point. This is easily realised by inspecting the condition ^5^^, = (for = 1). The solution 
of this condition is greatly simplified by the use of the Boltzmann approximation and the 
assumption that isospin symmetry leads to the approximate solution A„ = = A^ and 
lu — ld = Iq- Neglecting trivial solutions, where 7^ = 0, the zero strangeness condition can 
be solved to give 

_ FK{T)-FH{T)^q{\ + \-^) 

2F^{T) ■ ^^^> 

In eq. (22), Fx corresponds to the Kaon mesons, Fh to the Hyperon baryons (A's and 
E's) and Fs to the S baryons, while the summation 

Fm = ^,i:9^mlK,(^) (23) 

to all particles i of the same family is implied. In the above relation, K denotes a modified 
Bessel function of the second kind. It is evident from eq. (22) that the increase of the relative 
chemical equilibrium factor for light quarks, 7^ and the increase of the light quark fugacity, 
A^, leads, at constant temperature, to the decrease of factor 7^. 



4. Inclusion of pentaquarks 

There has been recent evidence of hadrons containing five quarks. These 5-quark states 
are the e+(1540) [22] with / = and quark content uudds and the 5*(1862) with / = 3/2. 
The content of the states S*(1862) is ssddu (for the state with electric charge Q=-2), ssudu 
(with Q=-l), ssudd (with Q=0) and ssuud (with Q=+l). The existence of the first three 
of the states S*(1862) has been confirmed [23]. The motivation to investigate the effect of 
the pentaquark states comes from the fact that these states can alter eqs. (l)-(6), since they 
introduce additional hadronic families, each of which is accompanied by completely different 
functions between the system fugacities compared to the ones in the already known families. 
This can be easily realised if the corresponding equation of (22) is written down as 

^ Fk{T) + Fe{T)^l{\l + A-^)(A, + A,'^) - Fh{T)^M, + K") .04^ 

2[Fs(T) + Fs*(T)72] " ^ > 

The existence of © hadron drives 7^ to higher values with a strong dependence on 7^ and 
Aq, whereas the inclusion of the S* states contributes to the decrease of 7^. 

The system of eqs. (l)-(6) is then solved with the inclusion of the ©"'"(1540) and 5*(1862) 
states with the same partition functions for the HG and the QGP phase and for the same 
parameters 5, Vq as in the case without the inclusion of the pentaquarks (lines (a) of the 
previous section). The resulting curves are lines (b), shown in Figs. 1-4. The adopted form 
of the ratio of volumes x again produces value X\ — 1.1 at A„ 1 = 14.2 (while e = 1.5) and is 
plotted in Fig. 5. 

Lines (a) and (b) record the difference induced in the transition curve by the inclusion 
of the pentaquarks, if parameters B and Vq remain the same and x acquires the same value 
at a given value of A^. However, none of these parameters is known and so the difference 
in the transition curve cannot provide evidence for the existence of pentaquarks. For this 
reason lines (c) are drawn in Figs. 1-4. These lines represent a solution for the transition 
curve without the inclusion of the pentaquarks but for a different choice of parameters {B 
remains the same, Vq — 1.29/(45), xi = 1.12 at A^ 1 — 14.2 and e = 1.15). It is evident, 
now, that lines (b) (with pentaquarks included) and lines (c) (without pentaquarks) almost 
coincide. 

5. Application to heavy- ion data 

The fact that the transition territory between the hadronic and the QGP phase is re- 
stricted to a line in the space of temperature and chemical potentials produces a direct 
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connection between the variables after the phase transition and the corresponding ones be- 
fore: they must coincide. This is not the case when the transition territory is aUowed to 
occupy a surface. Then the connection between the thermodynamic variables of the hadronic 
and the quark phase is broken. The system, as it crosses the transition territory to enter the 
state of hadrons, loses its "memory" of the plasma state. 

In a system where the hadronic state carries the memory of its preceding state, one may 
use the freeze-out variables as a diagnostic tool for a primordial QGP phase. Assuming that 
(a) a quark-gluon phase has been formed in a collision experiment and (b) the chemical 
freeze-out occurs right after the transition to the hadronic phase, then the freeze-out ther- 
modynamic variables have to fulfil constraints (l)-(6). If, on the contrary, no quark-gluon 
state is formed before hadronization, then the restriction on the freeze-out conditions of the 
system is diminished only to eqs. 



a set of constraints which will be referred to from now on as set A. 

The thermodynamic variables are extracted through a fit of the experimentally measured 
particle multiplicities or ratios to a statistical model. Such a technique has been successful. 

If now the additional constraints (l)-(6) are imposed, the question that arises is whether a 
successful fit is also produced or the restrictions that these constraints imply are inconsistent 
with the data. 

These ideas will now be used to analyse the freeze-out variables of four heavy-ion ex- 
periments. The application of eqs. (l)-(6) require the knowledge of the partition functions 
for the hadron and the quark phase. The particular functions used in sections 4 and 5 
to demonstrate a solution include the arbitrariness in the choice of the parameters B, Vq 
and quantity x. Since it is unwanted for the extracted variables to depend on the choice 
of unknown parameters, it is better to form and apply a subset of constraints which are 
completely independent from these parameters. 

When the system of eqs. (l)-(6) is valid, eqs. (3)-(6) can equivalently be rewritten as 



''^Bhg{T, {A, 7}) - 2(3 uq^^ 



(25) 



(26) 



^Shg {^^ 7}) - 2/3 Uq^^ 



(27) 



nB^^,(r,{A,7}) = 2P nQ^«,(r,{A,7}) 



(28) 



in 



nsHGiT,{^rf}) = , 



(29) 



^5«cp(r,{A,7}) = 0. 



(30) 



Eq. (27) results from cqs. (2), (3) and (5), whereas eq. (30) is due to eqs. (4) and (6). It is 
easily checked that the common denominator (see eq. (12)), which inchides the parameter 
Vq, cancels out from both sides of eq. (27). The same is true for eq. (29), while eq. (30) is 
easily solved to give = 1. Eqs. (27)- (30) (referred to as Set B) now form a set of equations 
that do not depend on the parameters Vq for the particle size, B (MIT bag constant) nor 
the ratio x applicable to the first order transition line. 

On the contrary, eqs. (1) and (2) are model dependent and contain unknown parameters. 
However, if the freeze-out parameters are determined, they can be inserted in eq. (2) to 
determine Vq (assuming that x is known) and then eq. (1) can be used to determine B. This 
task serves to show that eqs. (1) and (2) have a real solution and contributes to the overall 
consistency of the technique. 

In the following the freeze out variables for the experiments S + S [24] , S + Ag [25] 
(NA35) at beam energy 200 AGeV, Ph + Ph [26] (NA49) at beam energy 158 AGeV and 
Au -\- Au [27] (STAR) at -^snn — 130 GeV will be analysed. The data used are hsted in 
Table 1 and they are in all cases full phase space multiphcities, except from the RHIC data 
which are measured in the midrapidity. The experiments to be analysed are chosen because 
they do not produce great baryon chemical potential at freeze out and so they are probably 
at the crossover area [20], allowing one to set x = 1. The technique can be applied to the 
first order transition case, determining the freeze-out variables, since the equations of set B 
do not depend on x, but then the parameters Vq and B cannot not be uniquely determined. 

The theoretical calculation of the particle multiplicity necessary to perform a fit to the 
experimental data has been carried out with the partition function (7), (10). The right 
Bose-Einstein or Fermi-Dirac statistics for every particle has been used throughout the cal- 
culations. The feeds from the decay of resonances have also been included. The value of (3 
is set to 1 in the case oi S -\- S, 1.1 in the case oi S -\- Ag, 1.27 for Ph -\- Ph and 1.25 for 
Au + Au. 

In Tables 2 (a)-(b) the freeze-out variables are extracted with the assumption of no 
primordial QGP phase (constraints of set A), whereas in Tables 3 (a)-(b) the primordial 
QGP phase is assumed, thus applying constraints of set B. Another matter concerning the 
analysis is the observation that the inclusion of the pion multiplicity deteriorates the fit 
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[8,28]^. Since the quality of the fitting procedure is of importance in evaluating the results 
and a bad fit, when the constraints of set B are imposed, may be partly due to the presence 
of the pion multiplicity, two fits are performed in each case, one with all the multiplicies 
included (Tables (a)) and one without the multiplicities that contain pions (Tables (b)). 
This makes clear the effect of the pion multiplicities in the overall procedure. In any case, 
the fits without the pions are in general more reliable. 

The first observation which can be made by comparing the first and the second part 
of every table is that the inclusion of pentaquarks has negligible effect on the evaluated 
parameters or the quality of the fit. So one can safely draw equivalent conclusions by 
performing the analysis with the pentaquarks or without them. 

For the S + S and S + Ag data the fit with set B is of medium quality {x^/dof = 2.95 
and 1.92, respectively) when the pions are present. This is not far worse, though, than the 
fit in the case of set A. When the pions are excluded, the fit with set B turns out to be very 
good {x^/dof — 0.47 and 0.065, respectively), while the temperature remains at acceptable 
values, proving these cases to be completely compatible with a primordial quark-gluon phase. 
Another observation is that the imposition of set B with respect to set A, in every case, does 
not produce a dramatic change in the quality of the fit. 

In the case of Pb + Pb the fit is relatively good with the imposition of set A {x^/dof — 
2.50 — 1.79), but the imposition of set B severely worsens the quality of the fit. The situation 
cannot be remedied with the exclusion of pions and x^/dof remains at the value of 18.0—18.2. 
The conclusion to be drawn from the bad quality of the fit with the imposition of the 
constraints of set B, is that the data of this experiment are not compatible with a preceding 
QGP phase. The dramatic change in the quality of the fit between the cases of set A and B, 
should be noted as well. Also the freeze-out temperature in case of set B is unreahstically 
high and rises enormously with respect to the case of set A. 

The findings concerning the S + S and S + Ag data are also in agreement with the 
proximity of the chemical freeze out points of these experiments to the Statistical Bootstrap 
critical line which was found in [8] . On the contrary, the freeze out point of Pb + Pb was 
not found to possess such an attribute in [30], which is also in agreement with the present 
results. 

In case of RHIC and in case of set A the vahic of x^/dof is 1.51, when the pions are 
^The presence of excess of pions, though, can be connected with a primordial high entropy phase or with 
a phase with the chiral symmetry restored [29]. 
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included and this value diminishes to 0.229, when the pions are excluded. Similar results are 
obtained when set B is imposed. The fit in the presence of pions is of not so good quality 
{x^/dof — 3.86) and the temperature aquires too high value. The fit turns out to be quite 
good, though, when the pions are excluded {x^/dof — 1.19) and the temperature remains 
at acceptable values, so the thermodynamic parameters are compatible with a QGP phase. 

The extracted parameters in case of set B are inserted to eqs. (1) and (2) and the 
parameters Vq and B^^^ are also determined. It is interesting that in the cases oi S + S, 
S + Ag and Au + Au (without the pions) , which have been proven to be compatible with 
set B, all the calculated values of Vq and S^/^ are close, compatible with a unique value for 
these parameters. On the contrary, the Ph + Ph case produces values of Vq and B^/^ with 
no connection with the rest cases. 

The necessity of the expansion of the fugacity sector with the partial equilibrium fugacitics 
is also revealed with the application of the present technique. If these fugacities are set to 
lu — Id — Is — '^1 then the sector of the phase space which is compatible with the QGP- 
hadron transition is severely limited. In such a case, if a similar fit to the one with set B 
is performed, apart from the fact that eqs. (27) and (28) cannot be accommodated, the fit 
turns out to be worse. The result in the case without the pions is then jdof = 0.617, 1.09, 
28.0 and 1.83 ior S + S, S + Ag, Ph + Ph and Au + Au respectively. Then in the case of 
Au + Au the compatibility with the QGP phase turns out to be dubious. 

Another general observation, which can be drawn from Tables 2-3, is that in the cases 
which are compatible with a primordial QGP phase {S + S + Ag and Au + Au) the 
inclusion of the pions in the fits produces a dramatic increase in the /dof value relative to 
the fits without the pions. Also the fitted temperature is calculated to be much higher in the 
presence of pions. This is checked from comparing the value of jdoj and the temperature 
in Tables (b) and (a). Although in the present work the fugacities 7„ and 7^ are used, 
which describe off chemical equilobrium effects and the pion content is u and d quarks, this 
dramatic increase in the jdof and the temperature when the pions enter the fit, persists. 
Thus, for the QGP compatible cases, the 7„ and 7^ fugacities cannot improve the fit to 
acceptable limits and reveal that an excess of pions is present to these cases. 

It is evident by incorporating all the previous analyses, that for the S + S, S + Ag and 
Au + Au cases, the freeze-out thermodynamic variables should be considered these of Table 
3 (b). For the Pb + Ph case the relevant freeze-out variables should be considered these of 
Table 2 (a) or 2 (b). 



Gathering the observations made in this section, what it found is that the freeze-out 
conditions of some experiments are compatible with a QGP phase (constraints B). This is 
revealed from the fact that the imposition of constraints B produces fits with acceptable 
values of jdoj and temperature. The opposite is true for the experiments incompatible 
with the QGP phase. In such cases the imposition of constraints B leads to high values 
/dof and temperature, though these values are acceptable, when only constraints A (no 
QGP phase assumed) are imposed. Also, the QGP compatible cases present an excess in the 
pion multiplicity which cannot properly be fitted (in these cases) by the 7^ and fugacities. 
This is seen by a great increase in the value of jdoj and temperature when the pions are 
included in the fit in comparison to the fits without pions, for both cases of constraints A 
and B. 

6. Conclusions 

Although two different partition functions are used for the description of the quark and 
hadronic side of matter, it is possible to preserve the continuity of all chemical potentials and, 
of course, temperature and pressure (Gibbs equilibrium conditions) at the transition region, 
which is confined to a curve. Also, all the constraints imposed by the conservation laws of 
quantum quantities can be applied, leading, at the same time, to a non-trivial solution of the 
thermodynamic variables into a three quark flavour system. The key issue for the success of 
this project is the expansion of the fugacity sector of the available variables and the, already, 
introduced relative chemical equilibrium variables can be used to serve that purpose. 

Despite the fact that the space of the thermodynamic variables is extended, the restric- 
tions imposed on the transition points produce relations among these variables. A part of 
these relations, in a simple form, is expressed in eq. (22) or (24). 

The restrictions on the freeze-out conditions, imposed by the existence of a quark-gluon 
state in the early stages after a collision experiment, can be applied to every case where 
the thermalisation of the produced hadrons has been proven. They can serve to separate 
the experiments compatible with QGP state from those which are not. The expansion of 
the fugacity sector with the partial equilibrium fugacities, though, magnifies the part of the 
phase space allowed by such constraints. 

The whole methodology which was presented can be used for every grand canonical 
partition function adopted for the description of the HG or QGP phase. The inclusion of 
interaction is crucial for the prediction of the critical point and the volume expansion ratio. 
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which could not be determined by the models used in this work. At the moment, lattice 
calculations have led to the determination of the accurate quark-gluon equation of state 
with three quark flavours at finite baryon chemical potential [20,31]. It would be interesting, 
though, if these calculations could be extended with the inclusion of the relative chemical 
equilibrium variables for light and strange quarks, allowing for matching with the existing 
hadron gas models. For the hadronic side of matter the inclusion of the attractive part of 
interaction can be incorporated via the statistical bootstrap [7,8], where the prediction of a 
critical point is also possible [32] . The incorporation of the full set of parameters 7^ to these 
studies would allow for a more precise matching with a primordial quark phase. 

Acknowledgement I would like to thank N. G. Antoniou, C. N. Ktorides and F. K. Di- 
akonos for fruitful discussions. 
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Tables 



s + s 


S + Ag 


Pb + Pb 


Au + Au 


K+ 


12.5 ±0.4 




15.5 ± 1.5 


Np 


362 ±5.1 


A 


17.20 ± 1.75 


K- 


6.9 ±0.4 


A 


15.2 ± 1.2 


K+ 


103 ±7.1 


A 


12.15 ± 1.25 


Ks' 


10.5 ± 1.7 


A 


2.6 ±0.3 


K- 


51.9 ±3.6 




2.11 ±0.23 


A 


9.4 ± 1.0 


P 


2.0 ±0.8 




81 ±4 




1.77 ±0.19 


A 


2.2 ±0.4 


p — p 


43 ±3 


</> 


7.6 ±1.1 


VL+Vt 


0.585 ±0.150 


P 


1.15 ±0.40 


B-B 


105 ± 12 


A 


53 ±5 


P 


26.37 ±2.60 


p-p 


21.2 ± 1.3 




186 ± 11 


A 


4.64 ±0.32 


P 


18.72 ± 1.90 


B-B 


54 ±3 








4.45 ± 0.22 




36.7 ±5.5 




98 ±3 








0.83 ±0.04 


<t> 


5.73 ±0.78 










VL 


0.62 ±0.09 


K*^ 


10.0 ±2.70 










Vl 


0.20 ±0.03 


TT + W 


239 ± 10.6 












619 ±35.4 




239 ± 10.6 










7r-(*) 


639 ± 35.4 


K+/K- 


1.092 ± 0.023 
















0.92 ±0.27 














Vl/Vl 


0.95 ±0.16 



This multiplicity has not been used in the fits where the pions are excluded. 



Table 1. The full phase space multiplicities from the coUision experiments S + S (NA35), 
S + Ag (NA35) and Pb + Pb (NA49), as well as the midrapidity multiplicities and ratios 

from Au + Au (STAR), used in the fits. 
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a) set A, Fit with all 




s + s 


S + Ag 


Pb + Pb 


Au + Au 


No pentaquarks 






/I 1 8 /I 
41.10/ 1 


1 7 Ki /7 
1 / .01/ / 


1 Q «0 /Q 
lo.DZ/ y 


1 yivi t'V j 


904 9 -1- Q 7 


oqc _|_ 1 Q 


^Q'^ A-\- c\ 
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1. / o zn u. ±o 


i. I I zn u.oo 
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1 . / zn u.uy i 
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i.zu zn u.iz 




i.ooy m u.uuy 


1 ZL7 -1- n 1 Q 
1 ( ZL u. ly 


1 fiAfl -1- n 071 
i .u^u zn u.u ( i 


i.uoy zn u.ui4t 






1 090 -1- n nzLi 


1 1 fi7 -1- n ni 
1 . ID / zn u.uio 


1 nni zL -1- n c\c\A7 
i.uui^ zn u.uu^/ 




0.567 ±0.077 


0.46 ±0.14 


0.448 ± 0.059 


0.164 ±0.091 


Id 


1.06 ±0.14 


0.70 ±0.17 


0.583 ±0.053 


0.606 ± 0.084 


Is 


0.517 ±0.024 


0.338 ± 0.043 


0.378 ±0.019 


0.284 ±0.069 




173 ±43 


350 ± 150 


3671 ± 400 


550 ± 180 


piB (MeV) 


234 ± 29 


315 ± 80 


297 ±21 


84 ±38 


With pentaquarks 


x'/dof 


6.20/3 


4.11/1 


17.49/7 


13.62/9 


T (MeV) 


205 ± 10 


235 ± 19 


193.7 ±3.0 


326 ± 76 




1.76 ±0.15 


1.77 ±0.33 


1.724 ±0.097 


1.20 ±0.12 


Xd 


1.340 ±0.070 


1.46 ±0.19 


1.639 ±0.071 


1.039 ±0.014 




1.022 ±0.022 


1.022 ±0.040 


1.168 ±0.013 


1.0015 ± 0.0047 


lu 


0.563 ± 0.075 


0.46 ±0.14 


0.446 ± 0.059 


0.164 ±0.091 


Id 


1.05 ±0.14 


0.70 ±0.16 


0.582 ± 0.053 


0.606 ± 0.084 


Is 


0.514 ±0.017 


0.340 ± 0.042 


0.376 ±0.019 


0.284 ±0.069 




173 ± 43 


350 ± 150 


3670 ± 400 


550 ± 180 


piB (MeV) 


235 ± 30 


313 ± 79 


297 ±21 


84 ±38 



1Q 



(b) set A, Fit with no tt's 





S + S 


S + Ag 


Pb + Pb 


Au ± Au 


No pentaquarks 


X /"'Oj 


U.OD/ z 


n /n 
u/ u 


o.y I/O 


1.00/ 


1 \1V1 t-V j 




zoy.o 


± I ^.o zn ^.u 


±o±.\j zn 0.^ 
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1. / oi zt u.UoZ 


i.DDDD 


1 fiQ7 _|_ n ORO 
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i.U I OD zt U.UUOi 




i.oi^ m u.uuu 


i .o ioU 


1 71 n -1- n OfiQ 
i. ( iu ZL u.uuy 


1 079Q -1- 00^;^; 
i.U ( zy zn u.uuoo 






u.yui 1 


1 1 fiQ -1- n 01 9 
i.iuy zn u.uiz 


1 ni 79 -1- 009/1 
i.ui / z zn u.uuz^ 




0.548 ± 0.027 


0.5819 


0.670 ± 0.048 


1.37 ±0.14 


Id 


0.719 ±0.061 


0.5377 


0.759 ± 0.055 


1.75 (fixed) 


Is 


0.588 ±0.022 


0.4650 


0.575 ± 0.063 


1.93 ±0.17 


yrpi 


67.3 ±9.2 


158.08 


2426 ± 30 


275 ± 23 


liB{MeV) 


365 ± 41 


442.4 


279 ± 17 


32.4 ± 1.9 


With pentaquarks 


x'/dof 


0.35/2 


0/0 


9.10/5 


1.85/8 


T{MeV) 


258 ± 23 


265.1 


176.1 ±5.0 


151.8 ±3.4 




1.773 ±0.084 


1.6664 


1.700 ±0.062 


1.0761 ±0.0052 


Xd 


1.522 ± 0.080 


1.8210 


1.705 ± 0.069 


1.0732 ± 0.0055 




0.943 ±0.018 


0.9612 


1.170 ±0.012 


1.0178 ±0.0025 


lu 


0.557 ±0.071 


0.5676 


0.646 ± 0.055 


1.36 ±0.13 


Id 


0.72 ±0.10 


0.5237 


0.741 ± 0.054 


1.75 (fixed) 


Is 


0.597 ±0.056 


0.4544 


0.555 ±0.065 


1.92 ±0.17 




67± 14 


153.42 


2490 ± 85 


274 ± 23 


HB{MeV) 


365 ± 44 


453.1 


282 ± 18 


32.6 ± 1.9 



Table 2. The results of fits on the S + 5 (NA35), S + Ag (NA35), Pb + Pb (NA49) and 
Au + Au (STAR) data with the imposition of the set of constraints A, without and with 
the inclusion of the pentaquark states. In part (a) all the multiplicities have been included 
in the fit and in part b) the multiplicities that contain pions have been excluded from the 

fit. In part (b) the fit on the STAR data has been performed with the variable 7^ held 
fixed at the given value. Table 2 (a) or 2 (b) contains the extracted parameters which are 
concluded as the freeze-out conditions for the Pb + Pb case. 
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(a) set B, Fit with all 
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Pb + Pb 
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Is 


0.401 ± 0.094 


0.308 ± 0.095 


0.188 ±0.011 


0.322 ± 0.061 


yrpZ 


152.2 ±9.5 


280 ± 12 


678 ± 77 


378 ± 53 


/iB {MeV) 


313 ±36 
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702 ± 74 


84 ±22 


VQ{lQ-^^MeV-^) 


5.4 ±2.0 


4.9 ± 1.4 


2.33 ±0.35 


3.58 ±0.62 


B^l"^ {MeV) 


325 ± 37 
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548 ± 40 


441 ± 83 


With pentaquarks 


xVdof 


14.99/5 


5.78/3 


162.65/9 


42.53/11 


T {MeV) 


246 ± 22 


278 ± 74 


439 ± 32 


346 ± 63 




1.536 ±0.038 


1.613 ±0.046 


1.668 ±0.084 


1.082 ±0.021 


Xd 
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Id 
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0.322 ±0.061 




149.1 ± 8.5 


274 ± 11 
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HB {MeV) 
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85 ±22 
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5.5 ± 1.5 


4.9 ± 1.6 


2.31 ± 0.35 


3.57 ±0.63 


B^l"^ {MeV) 


328 ± 31 


362 ± 100 


551 ±41 


442 ± 82 
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(b) set B, Fit with no tt's 
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Is 


0.847 ±0.071 


0.601 ± 0.042 


0.199 ±0.012 


0.81 ± 0.27 


yrpZ 


94.1 ±6.5 


199 ± 10 


621 ± 67 


269 ± 59 


/iB {MeV) 


275.3 ±9.3 


327.0 ± 8.4 


778 ± 79 


49.5 ±8.1 


VQ{lQ-^^MeV-^) 


6.8 ± 1.0 


6.2 ±0.8 


2.2 ± 0.4 


5.4 ±4.6 


B^l"^ {MeV) 


278.9 ± 7.3 


290.7 ±7.6 


557 ± 48 


301 ± 41 


With pentaquarks 


xVdof 


1.55/4 


0.13/2 


128.04/7 


10.80/9 


T {MeV) 


196.7 ±4.6 


216.3 ±5.2 


446 ± 38 


221 ±26 




1.608 ±0.023 


1.663 ±0.011 


1.746 ±0.073 


1.075 ±0.012 


Xd 


1.603 ±0.023 


1.694 ±0.012 


1.817 ±0.083 


1.081 ±0.012 


7« 


0.948 ± 0.069 


0.726 ± 0.044 


0.261 ±0.018 


0.69 ±0.20 


7rf 


0.956 ±0.070 


0.746 ± 0.045 


0.280 ±0.020 


0.74 ±0.21 


7s 


0.865 ±0.078 


0.562 ±0.039 


0.198 ±0.012 


0.79 ±0.25 




88.5 ± 6.9 


200 ± 10 


615 ± 66 


272 ± 58 


HB {MeV) 


279.1 ±9.1 


338.2 ± 8.7 


781 ± 80 


50.4 ±8.2 


VQ{lQ-^^MeV-^) 


7.2 ±0.9 


6.3 ±0.7 


2.2 ± 0.4 


5.4 ±4.3 


B^l"^ {MeV) 


282.1 ± 7.5 


298.4 ± 7.8 


560 ± 48 


305 ± 40 



Table 3. The results of fits on the 5 ± 5 (NA35), S + Ag (NA35), Pb + Pb (NA49) and 
Au + Au (STAR) data with the imposition of the set of constraints B, without and with 
the inclusion of the pentaquark states. In part (a) all the multiplicities have been included 
in the fit and in part (b) the multiplicities that contain pions have been excluded from the 
fit. Table 3 (b) contains the extracted parameters which are concluded as the freeze-out 
conditions for the S + S, S + Ag and Au + Au cases. 
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Figure Captions 

Fig. 1 Temperature as a function of the baryon chemical potential at the QGP-Hadron gas 
transition hne for B^/^ = 280 MeV, without (lines (a),(c)) and with (hne (b)) the inclusion 
of the pentaquark states. Lines (a) and (b) are calculated with Vq = 1.4/(45) and line (c) 
with Vq — 1.29/(45). Line (d) is the phase transition curve calculated from lattice QCD in 
[21]. 

Fig. 2 Relative chemical equilibrium variable of zi-quark, 7„, as a function of the baryon 
chemical potential at the QGP-Hadron gas transition hne. Lines (a)-(c) correspond to hues 
(a)-(c) of Fig. 1. 

Fig. 3 Relative chemical equilibrium variable of d-quark, 7^;, as a function of relative chemi- 
cal equilibrium variable of w-quark, 7„, at the QGP-Hadron gas transition line for the isospin 
symmetric case. Lines (a)-(c) correspond to lines (a)-(c) of Fig. 1. The line 7^ = 7„ is also 
depicted. 

Fig. 4 Relative chemical equilibrium variable of s-quark, 7^, as a function of the baryon 
chemical potential at the QGP-Hadron gas transition hne. Lines (a)-(c) correspond to hues 
(a)-(c) of Fig. 1. 

Fig. 5 Volume expansion ratio x — Vhg/Vqgp between pure hadron and pure QGP phase 
at the same transition point, as a function of the baryon chemical potential, which was used 

in calculations in Figs. 1-4. Lines (a)-(c) correspond to lines (a)-(c) of Fig. 1. In lines (a) 
and (b) xi = 1.1 at A„i = 14.2 and e = 1.5 and in line (c) xi = 1.12 at A^i = 14.2 and 
e= 1.15. 
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